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Local and nonlocal contributions to the total strain rate tensor Sy at any point x in a flow are 
formulated from an expansion of the vorticity field in a local spherical neighborhood of radius R 
centered on x. The resulting exact expression allows the nonlocal (background) strain rate tensor 
S^(x) to be obtained from SV,(x). In turbulent flows, where the vorticity naturally concentrates into 
relatively compact structures, this allows the local alignment of vorticity with the most extensional 
principal axis of the background strain rate tensor to be evaluated. In the vicinity of any vortical 
structure, the required radius R and corresponding order n to which the expansion must be carried 
are determined by the viscous lengthscale A^. We demonstrate the convergence to the background 
strain rate field with increasing R and n for an equilibrium Burgers vortex, and show that this 
resolves the anomalous alignment of vorticity with the intermediate eigenvector of the total strain 
rate tensor. We then evaluate the background strain field Sfj (x) in DNS of homogeneous isotropic 
turbulence where, even for the limited R and n corresponding to the truncated series expansion, the 
results show an increase in the expected equilibrium alignment of vorticity with the most extensional 
principal axis of the background strain rate tensor. 

PACS numbers: 47.27.-i,47.32.C-,47.27.De 



I. INTRODUCTION 

Vortex stretching is the basic mechanism by which ki- 
netic energy is transfered from larger to smaller scales in 
three-dimensional turbulent flows IlL 

SHU- An under- 
standing of how vortical structures are stretched by the 
strain rate field Sij (x) is thus essential to any description 
of the energetics of such flows. Over the last two decades, 
direct numerical simulations (DNS) [j| H, 0] and experi- 
mental studies [t| [13, HH, E2] of the fine-scale structure 
of turbulence have revealed a preferred alignment of the 
vorticity with the intermediate eigenvector of the strain 
rate tensor. This result has been widely regarded as sur- 
prising. Indeed the individual components of the inviscid 
vorticity transport equation, in a Lagrangian frame that 
remains aligned with the eigenvectors of the strain rate 
tensor, are simply 

Dwi Dw 2 Dw 3 

where s±, S2 and S3 are the eigenvalues of Sij. For in- 
compressible flow, si + S2 + S3 = 0, and then denoting 
si > s 2 > S3 requires si > and s 3 < 0. As a con- 
sequence, (JIJ would predict alignment of the vorticity 
with the eigenvector corresponding to the most exten- 
sional principal strain rate si. Yet DNS and experimen- 
tal studies have clearly shown that the vorticity instead 
is aligned with the eigenvector corresponding to the in- 
termediate principal strain rate S2. 
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A key to understanding this result is that, owing to 
the competition between strain and diffusion, the vortic- 
ity in turbulent flows naturally forms into concentrated 
vortical structures. It has been noted, for example in 
Refs. @, Gl P, that the anomalous alignment of the 
vorticity with the strain rate tensor Sij (x) might be ex- 
plained by separating the local self-induced strain rate 
field created by the vortical structures themselves from 
the background strain field in which these structures re- 
side. The total strain rate tensor is thus split into 

^■(x)=S«(x) + S?(x), (2) 

where S^ is the local strain rate induced by a vortical 
structure in its neighboring vicinity, and Sj^ is the non- 
local background strain rate induced in the vicinity of 
the structure by all the remaining vorticity. The vortical 
structure would then be expected to align with the princi- 
pal axis corresponding to the most extensional eigenvalue 
of the background strain rate tensor Sfj (x) . 

In the following, we extend this idea and suggest a 
systematic expansion of the total strain rate field Sij (x) 
that allows the background strain rate field S^(x) to 
be extracted. Our approach is based on an expansion 
of the vorticity over a local spherical region of radius R 
centered at any point x. This leads to an exact opera- 
tor that provides direct access to the background strain 
rate field. The operator is tested for the case of a Burg- 
ers vortex, where it is shown that the local self-induced 
strain field produced by the vortex can be successfully re- 
moved, and the underlying background strain field can be 
increasingly recovered as higher order terms are retained 
in the expansion. The anomalous alignment of the vor- 
ticity with respect to the eigenvectors of the total strain 
field is shown in that case to follow from a local switch- 
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From (0 and ([6]), Sy(x) can be expressed [I5[ as an 
integral over the vorticity field as 




FIG. 1: Decomposition of the vorticity field in the vicinity 
of any point x into local and nonlocal parts; the Biot-Savart 
integral in (|8} over each part gives the local and nonlocal 
(background) contributions to the total strain rate tensor Sij 
at x. 



ing of the principal strain axes when the vortex becomes 
sufficiently strong relative to the background strain. Fi- 
nally, the operator is applied to obtain initial insights 
into the background strain S^(x) in DNS of homoge- 
neous isotropic turbulence, and used to compare the vor- 
ticity alignment with the eigenvectors of the total strain 
field and of this background strain field. 
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Sij(x) = — (eikirj + nejki) -^uJi{x.')d 3 x' 
87T J A r° 



(7) 



As shown in Fig.[Tl the total strain rate in ([7]) is separated 
into the local contribution induced by the vorticity within 
a spherical region of radius R centered on the point x 
and the remaining nonlocal (background) contribution 
induced by all the vorticity outside this spherical region. 
The strain rate tensor in ([7]) thus becomes 



r<R 



SF 



I \—]d 

Jr>R 



3 x',(8) 



where [• • • ] denotes the integrand in ([7]) . The nonlocal 
background strain tensor at x is then 



S5(x) = ^(x)-55(x) 



(9) 



The total strain tensor Sij (x) in © is readily evaluated 
via ([6]) from derivatives of the velocity field at point x. 
Thus all that is required to obtain the background (non- 
local) strain rate tensor S^{x) via © is an evaluation 
of the local strain integral (x) in (J5J) produced by the 
vorticity field a^(x') within r < R in Fig. Q] 



A. Evaluating the Background Strain Rate Tensor 



II. THE BACKGROUND STRAIN FIELD 

The velocity u at any point x induced by the vorticity 
field <*>(x) is given by the Biot-Savart integral 
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u(x) = -!- / w(x') x X * cPx' 
4tt/ a jx-x'^ 



(3) 



where the integration domain A is taken to be infinite or 
periodic. In index notation ([3]) becomes 

m(x) = i / ^(xO ^'l^ V , (4) 



■in 



where em- is the cyclic permutation tensor. The deriva- 
tive with respect to Xj gives the velocity gradient tensor 



-^-Ui(x) = -!- / e ak LOi(-K.') 
dXj Att J a 



5 M _ 2 rkVj 



d 3 x' 



(5) 



where r = |x — x'| and r m = x m — x' rn . The strain 
rate tensor Sij at x is the symmetric part of the velocity 
gradient, namely 



(6) 



The vorticity field within the sphere of radius R can 
be represented by its Taylor expansion about the center 
point x as 



w;(x')| 



r<R 



W;(x) + (x' m - X rn ) 



dbJl 



dx r , 



(10) 
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Recalling that x m - x' m = r m and using a h b lm , ci mn , . . . 
to abbreviate the vorticity and its derivatives at x, we 
can write (fTU|) as 



W '( x ')lr<_R = a l ~ Trahra + \^ m T n Cl T , 



(11) 



Substituting (If I p in the integral in © and changing 
the integration variable to r = x — x', the strain tensor 
at x produced by the vorticity in R is 



(x) 



3 f ( i \ Tk 



(12) 



ai - r m bi m + 



d v . 



This integral can be solved in spherical coordinates cen- 
tered on x, with n = rsm0cos4>, ri = r sin 9 sin</>, and 
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r 3 = rcos6 for r E [0,R], 9 £ [0,tt], and e [0,2tt). To 
integrate (fT2")l note that 



B. Practical Implementation 



r<R 



TkTj , 3 4?T 

^^crr = , 



(13a) 
(13b) 



r<R 



f r k r j r m r n ,3 27T 2 / r c /1 o \ 

/ — ^ d r = — i? {S mn Ojk+ (13c) 

8mj&kn ~t~ ^mk^jn) ■ 

The resulting local strain rate tensor at x is then 
i? 2 

^ij (p^y — Qmn {^ijl^mn ~\~ £jil&mn ^inl&mj (14) 

+ ) + 0(i? 4 ) , 

where the contribution from the en term in (|12[) is zero 
since e^; = ~ e ju- F° r the same reason the first two terms 
in also cancel, giving 

i? 2 



S/(x) 
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+ S ni ) + 0(i? 4 ) . 



Recalling that q„ m = q ?( 



tracting with the 8 and e in (|15p. gives 



d 2 u>i/dx m dx n , and con- 



«5W - 1 
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dxi \ 3m dx m 



(16) 



0(R 4 



(17) 



Note that e im i duii/dx m = (V x u>) . and 

V x oj = V x (V x u) = V (V • u) - V 2 u . 

so for an incompressible flow (V ■ u = 0) the local strain 
rate tensor at x becomes 

(>/l,y 0(i? 4 ). (18) 



From with S'y from (jTHJ) we obtain the background 



strain tensor as 
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S£(x) = S«(x) + ^V%(x) + 0(i? 4 ) . (19) 



The remaining terms in (|19p result from the higher-order 
terms in (jlip . The contributions from each of these can 
be evaluated in an analogous manner, giving 

i? 4 



&ij ( x ) — 



R 2 . 



+ 
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-V 2 V 2 



(20) 



5 i: ,-(x) 



(2n - 2)!(4n 2 - 1) v ; 

where the terms shown in (f2"0"|) correspond ton = 1,2,.... 
The final result in (|2"0"1) is an operator that extracts the 
nonlocal background strain rate tensor S-j at any point 
x from the total strain rate tensor Sij. For the Taylor 
expansion in (|10|) . this operator involves Laplacians of 
the total strain rate field Sij{x). 



When using (|20|) to examine the local alignment of any 
concentrated vortical structure with the principal axes of 
the background strain rate field (x) in which it resides, 
the radius R must be taken sufficiently large that the 
spherical region |x' — x| < R encloses essentially all the 
vorticity associated with the structure, so that its local 
induced strain rate field is fully accounted for. Gener- 
ally, as R increases it is necessary in (|2T?|) to retain terms 
of increasingly higher order n to maintain a sufficient 
representation of w(x') over the spherical region. Thus 
for any vortical structure having a characteristic gradi- 
ent lengthscale X v , it can be expected that R must be of 
the order of A„, and n will then need to be sufficiently 
large to adequately represent the vorticity field within 
this sphere. However, since the local gradient lengthscale 
in the vorticity field in a turbulent flow is determined by 
an equilibrium between strain and diffusion, the vortic- 
ity field over the lengthscale A„ will be relatively smooth, 
and thus relatively low values of n may suffice to give a 
usable representation of u;(x'). This is examined in the 
following Section. 



III. TEST CASE: BURGERS VORTEX 

The equilibrium Burgers vortex [l], [jj [l6| is formed 
from vorticity in a fluid with viscosity v by a spatially 
uniform, irrotational, axisymmetric background strain 
rate field Sf^ that has a single extensional principal strain 
rate S zz directed along the z axis, as shown in Fig. [2] 
This simple flow, often regarded as an idealized model 
of the most concentrated vortical structures in turbulent 
flows, provides a test case for the result in (f20|) . The 



combined strain rate field SV/(x) produced by the vor- 
tex and the background strain flow should, when applied 




FIG. 2: Equilibrium Burgers vortex with circulation V and 
strain-limited viscous diffusion lengthscale A„ in a uniform, 
irrotational, axisymmetric background strain rate field Sfj (x) . 
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in (|20]) . produce the underlying background strain field 
(S^S^SfJ = (-±,-±,l)S xz at all x when R -> oo 
and all orders n are retained. For finite R and n, the 
resulting SS (x) will reflect the convergence properties of 



S zz , and the corresponding principal strain axis will point 
in the z direction. The vorticity is then aligned with the 
most extensional eigenvector of Sij. This remains the 
case until the vortex becomes sufficiently strong relative 
to the background strain rate that s > S zz , namely 



A. Strain Rate Tensor 

The equilibrium Burgers vortex aligned with the exten- 
sional principal axis of the background strain rate field 
has a vorticity field 



a r 

w(x) = w x (r)z = --jexp I 

7T A„ 



-ar) I z , 



(21) 



where T is the circulation, A„ is the viscous lengthscale 
that characterizes the diameter of the vortex, 77 = rj\ v 
is the radial similarity coordinate, and the constant a 
reflects the chosen definition of \ v . Following fs], A„ is 
taken as the full width of the vortical structure at which 
uj z has decreased to one-fifth of its peak value, for which 
a = 4 In 5. When diffusion of the vorticity is in equilib- 
rium [9J with the background strain, then 



A„ = V8a ( — 

"22 



1/2 



(22) 



The combined velocity field u(x) produced by the vortex 
and the irrotational background strain is given by the 
cylindrical components 



u r (r, 9, z) 

u$(r, 9, z) 
u z (r,9,z) 



S zz 

r 1 



2n\ v 77 

S zz z . 



(23a) 

[l-exp(-my 2 )] , (23b) 
(23c) 



The combined strain rate tensor for such a Burgers 
vortex is thus 





" S zz /2 


QV 

D r6 







Sy( x ) = 




S zz /2 





(24) 










S zz 





where S%g is the shear strain rate induced by the vortex, 
given by 



r 

Sr6 ( x ) = To" 
7TAf 



2N 1 



(25) 



From (|24[) , S'y(x) has one extensional principal strain 
rate equal to S zz along the z axis, with the remaining 
two principal strain axes lying in the r-9 plane and cor- 
responding to the principal strain rates 



tS zz i \S^g\ 



(26) 



As long as the largest s in (|2"6")l is smaller than S zz , the 
most extensional principal strain rate si of Sij will be 



S r g\ > n S Z z j 



which from (1251) occurs wherever 



1 \ / 2n 1 3tt (T/\l\ 1 



(27) 



(28) 



At any 77 for which (|28|) is satisfied, the most extensional 
principal axis of the combined strain rate tensor Sij (x) 
will switch from the z direction to instead lie in the r- 
9 plane. Since the vorticity vector everywhere points in 
the z direction, wherever (|2"5|) is satisfied the principal 
axis of Sij that is aligned with the vorticity will switch 
from the most extensional eigenvector to the intermediate 
eigenvector. This alignment switching is purely a result 
of the induced strain field Sy (x) locally dominating the 
background strain field S?(x). 

The dimensionless vortex strength parameter 



n = 



^max 
a S zz 



(29) 



on the right-hand side of (|28|) characterizes the relative 
strength of the background strain and the induced strain 
from the vortical structure, where u) m ax is obtained from 
([H]) at T) = 0. For 



n < n* « 2.45. 



(30) 



the background strain rate S zz is everywhere larger than 
the largest s in (|26p . and thus no alignment switching 
occurs at any 77. For f2 > Q* , alignment switching will 
occur over the limited range of rj values that satisfy (|28[) . 
With increasing values of fl, more of the vorticity field 
will be aligned with the intermediate principal axis of 
the combined strain rate tensor, even though all of the 
vorticity field remains aligned with the most extensional 
principal axis of the background strain rate tensor. 

Figure [3] shows the vorticity uj z and the induced shear 
strain component — S% as a function of 77. The horizontal 
dashed lines correspond to three different values of f2, 
and indicate the range of rj values where the alignment 
switching in (j28|) occurs for each f2. Wherever — S^g is 
above the dashed line for a given f2, the vorticity will be 
aligned with the local intermediate principal axis of the 
combined strain rate field. 

In principle, regardless of the vortex strength parame- 
ter fi, at any 77 the result in (|2"0|) can reveal the alignment 
of the vorticity with the most extensional principal axis 
of the background strain field Sfj . However, this requires 
R to be sufficiently large that a sphere with diameter 2R, 
centered at the largest 77 for which —S" g in Fig. [3] is still 
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FIG. 3: Similarity profiles of io z (rf) and S r e(il) for any equi- 
librium Burgers vortex; wherever —S r e exceeds the horizontal 
line determined by the relative vortex strength parameter fl 
iu (|29p the most extensional principal axis of the total strain 
rate ^(x) switches from the z-axis to lie in the r-8 plane. 



above the horizontal dashed line, will enclose essentially 
all of the vorticity associated with the vortical structure. 
As ft increases, the required R will increase accordingly 
as dictated by (|2"5)) , and as R is increased the required n 
in (|2"0)) also increases. 

Irrespective of the value of fi, when (l2Tj"l) is applied to 



the combined strain rate field SV,(x) in (|2~4")) and 
if R = (i?/A„) — > oo and all orders n are retained then 
the resulting Sfj (x) should recover the background strain 
field, namely 



the approximation becomes increasingly poorer, but the 
1/r 2 decrease in the Biot-Savart kernel in neverthe- 
less renders it adequate to account for most of the vortex- 
induced strain rate field. At the largest radial location, 
corresponding to the bottom right panel of Fig. [4] the ap- 
proximation becomes relatively poor, however at large i] 
values the vortex-induced strain is sufficiently small that 
it is unlikely to lead to alignment switching for typical fl 
values. 

Figure [5] shows the shear component S^ e (rj) of the 
background strain rate tensor obtained via ([20]) for var- 
ious n and R as a function of r\. In each panel, the 
black curve shows the total strain rate (rj) and the col- 
ored curves show the background strain rate (r/) from 

(PCJ|) for the (n, R) combinations listed. The horizontal 
dashed line corresponding to ft = (3/2)17* reflects the rel- 
ative vortex strength, and shows the range of rj where the 
anomalous alignment switching occurs due to the vortex- 
induced strain field. Wherever the — S^g curves are above 
this line, the vorticity there will be aligned with the inter- 
mediate principal axis of the combined strain rate tensor 
Sij ■ Figure [SJa) examines the effect of increasing the 
radius R of the spherical region for fixed order n = 6. 
It is apparent that with increasing R the resulting —S^e 
converges toward the correct background strain field in 
(f3"Tj) . For the value of ft shown, it can be seen that for 
R > 0.5AI, the resulting S^ g is everywhere below the hor- 
izontal dashed line, indicating that the vorticity every- 
where is aligned with the most extensional principal axis 
of the resulting background strain rate tensor S*j (x) from 







(31) 



for all x, and the vorticity should show alignment with 
the most extensional principal axis of . For finite R/X v 
and various orders n, the convergence of from (|2"0|) to 
this background strain field is examined below. 



B. Convergence of the Background Strain 

The accuracy with which ([2"P|) can recover the back- 
ground strain field S§(x) that acts on a concentrated 
vortical structure depends on how well the expansion in 
(flT))) represents the vorticity field within the local spheri- 
cal neighborhood R. Figure 2] shows the results of a local 
sixth-order Taylor series approximation for the vorticity 
in (j2Tj) at various radial locations across the Burgers vor- 
tex. In each panel, the blue square marks the location x 
at which the sphere is centered, and the red dashed curve 
shows the resulting Taylor series approximation for the 
vorticity. On the axis of the vortex, the approximated 
vorticity field correctly accounts for most of the circula- 
tion in the vortex, and thus the induced strain field from 
the vortex will be reasonably approximated. Off the axis, 
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FIG. 4: Accuracy of the Taylor expansion for the local vortic- 
ity in (|10p for a Burgers vortex, showing results for 6th order 
approximation. In each panel, solid black curve shows actual 
vorticity profile, and red dashed curve gives approximated 
vorticity from derivatives at location marked by square. 
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FIG. 5: Convergence of background strain rate field Sfj(x.) for 
a Burgers vortex, obtained from total strain rate field Sij(x) 
using (|20[) for various (n,R) combinations, where R = R/X u . 
Shown are effects of increasing R for fixed n = 6 (top), in- 
creasing n for fixed R — 0.65 (middle), and increasing n and 
R simultaneously (bottom). The dashed horizontal lines fol- 
low from d27| and 



In Fig. \5j(b) similar results are shown for the effect of 
increasing the order n of the expansion for the vorticity 
field for fixed R = 0.65. It is apparent that the effect of 
n is somewhat smaller than for R in Fig.^et). Moreover, 
the results suggest that the series in (fSTTj) alternates with 
increasing order n. For this f2 and R, even n = 3 is 
seen to be sufficient to remove most of the vortex-induced 
shear strain, and thus reduce (x) below the horizontal 
dashed line. For these parameters, the S? 6 field from 
(|2T))) would thus reveal alignment of the vorticity with the 



most extensional principal axis of the background strain 
tensor throughout the entire field. 

Figure [5jc) shows the combined effects of increasing 
both R and n, in accordance with the expectation that 
larger R should require a higher order n to adequately 
represent the vorticity field within the spherical region. 
The shear strain rate field shows convergence to the cor- 
rect background strain field in pTj) . The convergence of 
the shear strain rate S^gfe) to zero in the vicinity of the 
vortex core is of particular importance. The systematic 
reduction in the peak remaining shear stress indicates 
that, even for increasingly stronger vortices or increas- 
ingly weaker background strain fields as measured by f2, 
the resulting S^ g (x.) from (j2"0)) will reveal the alignment 
of all the vorticity in such a structure with the most ex- 
tensional principal strain axis of the background strain 
field. 



IV. VORTICITY ALIGNMENT IN TURBULENT 
FLOWS 



Having seen in the previous Section how (f2T))) is able 
to reveal the expected alignment of vortical structures 
with the most extensional eigenvector of the background 
strain rate in which they reside, in this Section we apply 
it to obtain insights into the vorticity alignment in tur- 
bulent flows. In particular, we examine the alignment at 
every point x of the vorticity uj relative to the eigenvec- 
tors of the total strain rate tensor field Sij (x) and those 
of the background strain field data S?? (x). This analysis 
uses data from a highly-resolved, three-dimensional, di- 
rect numerical simulation (DNS) of statistically station- 
ary, forced, homogeneous, isotropic turbulence [l7l lis|. 
The simulations correspond to a periodic cube with sides 
of length of 2ir resolved by 2048 3 grid points. The Taylor 
microscale Reynolds number R\ is 107. 

The DNS data were generated by a pseudospectral 
method with a spectral resolution that exceeds the stan- 
dard value by a factor of eight. As a result, the highest 
wavenumber corresponds to k max r]K = 10, and the Kol- 
mogorov lengthscale r\K — v z ^ 4 / (e) 1 ^ 4 is resolved with 
three grid spacings. This superfine resolution makes it 
possible to apply the result in (|20f for relatively high or- 
ders n, which require accurate evaluation of high-order 
derivatives of the DNS data. In Schumacher et al. [TtJ 
it was demonstrated that derivatives up to order six are 
statistically converged. More details on the numerical 
simulations are given in Refs. (l7l. Il8j. 

Figure[S]gives a representative sample of the DNS data, 
where the instantaneous shear component S12 of the total 
strain rate tensor field Sij (x) is shown in a typical two- 
dimensional intersection through the 2048 3 cube. The 
data can be seen to span nearly 700 t]k in each direction. 
The 512 2 box at the lower left of Fig. [6] is used here to 
obtain initial results for alignment of the vorticity with 
the eigenvectors of the background strain rate tensor. 

The background strain rate tensor field (x) is first 
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FIG. 6: Instantaneous snapshot of total strain rate com- 
ponent field 5i2(x) in a two-dimensional slice through a 
highly-resolved three-dimensional (2048 3 ) DNS of homoge- 
neous, isotropic turbulence [rH [Isl]. Axes are given both in 
grid coordinates (i — 1 . . . 2048) and normalized by the Kol- 
mogorov length wk ■ Box indicates region in which background 
strain rate field Sfj (x) is computed in Fig. [7] 



extracted via (|2"tJ)) from SV,(x) for n — 3 and vari- 
ous (R/t)k)- Higher-order evaluation of the background 
strain rate is not feasible, as the results in Ref. [l7( show 
that only spatial derivatives of the velocity field up to 
order six can be accurately obtained from these high- 
resolution DNS data. For n = 4, the expansion in (fSTTj) in- 
volves seventh-order derivatives of the velocity field, and 
the background strain evaluation becomes limited due to 
the grid resolution. The results are shown and compared 
in Fig. where the shear component S±2 of the full strain 
rate tensor is shown at the top, and the corresponding 
nonlocal (background) component Sf 2 and local compo- 
nent S^2 are shown, respectively, in the left and right 
columns for (R/r/x) = 2.5 (top row), 3.5 (middle row), 
and 4.5 (bottom row). Consistent with the results from 
the Burgers vortex in Fig. [5J as (R/rix) increases the 
magnitude of the extracted local strain rate in the right 
column increases. However, for the largest (R/r/x) = 4.5 
case, n = 3 appears to be too small to adequately rep- 
resent the local vorticity field. This leads to truncation 
errors which are manifested as strong ripples in the back- 
ground and local strain fields (see panels (/) and (g))- 

The results in Fig. [7] thus indicate that radii up to 
(R/t)k) = 3.5 in combination with n = 3 can be used to 
assess alignment of the vorticity vector with the eigenvec- 
tors of the background strain rate field. Figure [5] shows 
the probability densities of the alignment cosines for the 
vorticity vector with the total strain rate tensor and with 
the background strain rate tensors from PU|). We com- 
pare Sij (Fig. |Hb) with for (R/m) = 2.5, n = 3 (Fig. 
W) and for (R/n K ) = 3.5, n = 3 (Fig.|5]c). The re- 
sults for alignment with the total strain rate tensor are 
essentially identical to the anomalous alignment seen in 



numerous other DNS studies 0, H, 0] and experimental 
studies [1, H, [l(| EH EH , which show the vorticity to be 
predominantly aligned with the eigenvector correspond- 
ing to the intermediate principal strain rate. However, 
the results for the Burgers vortex in the previous section 
show that such anomalous alignment with the eigenvec- 
tors of the total strain rate tensor is expected when the 
local vortex strength parameter f2 is sufficiently large to 
cause alignment switching. 

By comparison, the results in Fig. [5] (b) and (c) ob- 
tained for the alignment cosines of the vorticity vector 
with the background strain rate tensor Sj^ from (|20|) 
show a significant decrease in alignment with the inter- 
mediate eigenvector, and an increase in alignment with 
the most extensional eigenvector. While data in panel 
(6) show only a slight change compared to those in (a), 
the results in panel (c) demonstrate that our decomposi- 
tion can indeed diminish the anomalous alignment signif- 




FIG. 7: Total strain rate component field 5i2(x) (a), with cor- 
responding results from (|20[) for nonlocal (background) field 
5f 2 (x) (left) and local field Sf 2 (x) (right) for (R/vk) = 2.5 
(b,c), (R/vk) = 3.5 (d,e), and (R/m) = 4.5 (f,g), all with 
n = 3. 
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FIG. 8: Probability densities of alignment cosines for the vorticity with the eigenvectors of the strain rate tensor, showing 
results for Sij (a) and for Sfj using (R/t]k) = 2.5 with n = 3 (b) and (R/t]k) = 3.5 with n = 3 (c). 



icantly. This is consistent with the results for the Burgers 
vortex in the previous Section, and with the hypothesis 
that the alignment switching mechanism due to the local 
contribution SR to the total strain rate tensor is the pri- 
mary reason for the anomalous alignment seen in earlier 
studies. It is also consistent with the expected equilib- 
rium alignment from ([1]). While a more detailed study is 
needed to examine possible nonequilibrium contributions 
to the alignment distributions associated with eigenvec- 
tor rotations of the background strain field, as well as 
to definitively determine the R and n convergence of the 
background strain rate tensor in Fig. the present find- 
ings support both the validity of the result in ([20]) for 
extracting the background strain rate tensor field <Sy (x) 
from the total strain rate tensor field Sij(x), and the hy- 
pothesis that at least much of the anomalous alignment 
of vorticity in turbulent flows is due to the differences be- 
tween the total and background strain rate tensors and 
the resulting alignment switching noted herein. 



V. CONCLUDING REMARKS 

We have developed a systematic and exact result in 
([20| that allows the local and nonlocal (background) con- 
tributions to the total strain rate tensor Sij at any point 
x in a flow to be disentangled. The approach is based on 
a series expansion of the vorticity field in a local spherical 
neighborhood of radius R centered at the point x. This 
allows the background strain rate tensor field (x) to 
be determined via a series of increasingly higher-order 
Laplacians applied to the total strain rate tensor field 
Sij (it). For the Burgers vortex, with increasing radius 
R relative to the local gradient lengthscale A„ and with 
increasing order n, we demonstrated convergence of the 
resulting background strain tensor field to its theoreti- 
cal form. We also showed that even with limited R and 
n values, the local contribution to the total strain rate 
tensor field can be sufficiently removed to eliminate the 



anomalous alignment switching throughout the flow field. 
This conclusion is expected to also apply to the more re- 
alistic case of a non-uniformly stretched vortex where 

s„=f{z) Earn nap. 

Consistent with the results for the Burgers vortex, 
when (|20p was used to determine the background strain 
rate tensor field (x) in highly-resolved DNS data for a 
turbulent flow, the anomalous alignment seen in previous 
DNS and experimental studies was substantially reduced. 
We conclude that (|20p allows the local background strain 
rate tensor to be determined in any flow. Furthermore, 
we postulate that the vorticity vector field in turbulent 
flows will show a substantially preferred alignment with 
the most extensional principal axis of the background 
strain rate field, and that at least much of the anoma- 
lous alignment found in previous studies is simply a re- 
flection of the alignment switching mechanism analyzed 
in Section III and conjectured by numerous previous in- 
vestigators. 

Lastly, the result in (j2"0|) is based on a Taylor series ex- 
pansion of the vorticity within a spherical neighborhood 
of radius R around any point x. Such an expansion inher- 
ently involves derivatives of the total strain rate tensor 
field, which can lead to potential numerical limitations. If 
larger R and correspondingly higher orders n are needed 
to obtain accurate evaluations of background strain rate 
fields, then otherwise identical approaches based on al- 
ternative expansions may be numerically advantageous. 
For instance, an expansion in terms of orthonormal basis 
functions allows the coefficients to be expressed as inte- 
grals over the vorticity field within r < i?, rather than as 
derivatives evaluated at the center point x. (For exam- 
ple, wavelets have been used to test alignment between 
the strain rate eigenvectors and the vorticity gradient in 
two-dimensional turbulence [Hi].) This would allow a re- 
sult analogous to ([20]) that can be carried to higher orders 
with less sensitivity to discretization error. The key con- 
clusion, however, of the present study is that it is possible 
to evaluate the background strain tensor following the 
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general procedure developed herein, and that when such 
methods are applied to assess the background strain rate 
fields in turbulent flows they reveal a substantial increase 
in the expected alignment of the vorticity vector with the 
most extensional principal axis of the background strain 
rate field. 
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